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Width of the QCD transition in a Polyakov-loop DSE model 
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We consider the pseudocritical temperatures for the chiral and deconfinement transitions within 
a Polyakov-loop Dyson-Schwinger equation approach which employs a nonlocal rank-2 separable 
model for the effective gluon propagator. These pseudocritical temperatures differ by a factor of 
two when the quark and gluon sectors are considered separately, but get synchronized and become 
coincident when their coupling is switched on. The coupling of the Polyakov-loop to the chiral 
quark dynamics narrows the temperature region of the QCD transition in which chiral symmetry 
and deconfinement is established. We investigate the effect of rescaling the parameter To in the 
Polyakov-loop potential on the QCD transition for both the logarithmic and polynomial forms of 
the potential. While the critical temperatures vary in a similar way, the width of the transition 
is stronger affected for the logarithmic potential. For this potential the character of the transition 
changes from crossover to a first order one when To < 210 MeV, but it remains crossover in the 
whole range of relevant To values for the polynomial form. 

PACS numbers: 11.10.Wx,12.38.Aw,12.38.Mh,12.39.Fe 



I. INTRODUCTION 

The QCD phase transition between highly excited 
hadronic matter and the quark-gluon plasma is presently 
under experimental investigation at ultra-relativistic 
heavy-ion collider facilities like the Relativistic Heavy- Ion 
Collider (RHIC) at the Brookhaven National Laboratory 
or the Large Hadron Collider (LHC) at CERN Geneva. 
Its theoretical description requires methods to solve QCD 
at finite temperature in the highly nonperturbative low- 
energy domain. At present, the only method to obtain 
ab-initio solutions of QCD in this domain is lattice QCD 
(LQCD). 

LQCD calculations are becoming available in the re- 
gion of physical quark masses [IHll and allow for a quanti- 
tative description of the equation of state (EoS) and a de- 
termination of the pseudocritical temperature. Although 
still afflicted with some uncertainty, those results can be 
used in phenomenological studies of the QCD transition 
relevant for current and future heavy ion experiments at 
RHIC and LHC. Hadron gas model calculations give a 
good description of the EoS in the confined phase. Fur- 
thermore, statistical models of hadron production [l(| 
can be used to extract the chemical freeze-out tempera- 
ture which serves as a lower limit for the deconfinement 
temperature. 

Even so, there are limitations to explore the full QCD 
phase diagram with LQCD. Hence dynamical models for 
the phase structure of QCD which can be calibrated with 
hadron phenomenology and with finite- T LQCD remain 
an important tool. Particularly useful are chiral quark 
models of the Nambu-Jona-Lasinio (NJL)- type [11H16} . 
but they suffer from nonphysical quark excitations at 
low temperatures, below T c . Also the contribution from 



the gluon sector to the thermodynamics is missing. A 
rather successful generalization has recently been sug- 
gested which fixes both these problems by coupling the 
chiral dynamics of the quark sector as modeled within the 
NJL model to a mean-field description of the gluon sector 
with an effective potential as a function of the Polyakov- 
loop (PL) variable, fitted to the pure gauge lattice sim- 
ulations for the Yang-Mills pressure [17fl . Within such a 
PNJL model []~sl - 420j , a remarkable agreement with other 
LQCD results as, e.g., for the chiral susceptibilities, could 
be achieved once the temperature in the fit of the gluon 
mean-field is appropriately scaled to the critical tempera- 
ture obtained in the lattice simulations [13, HH • The 
success of this type of chiral quark model in reproducing 
features of lattice QCD simulations has led to a number 
of applications as well as to extensions of the model. For 
example, we would like to mention here the extensions 
to include eight-quark interactions [23j and to consider 
a coupling (entanglement) between the scalar coupling 
constant and the Polyakov-loop [24[ . 

The deficiency of such a PNJL model is in the poor 
quark dynamics with, e.g., a constant, momentum inde- 
pendent quark mass function. Although the excitation 
of quark degrees of freedom is strongly suppressed at 
low temperatures by the destructive interference due the 
PL phase factors, a dynamical confinement mechanism 
is absent. In order to improve the quark dynamics of a 
PNJL model, the local current-current form of the quark 
interaction has been generalized to a nonlocal one [25l - 
[27j | . In this way, a dynamical quark mass function can 
be modeled using appropriate covariant ansatze for the 
form factor of the nonlocality. However, before a quan- 
titative comparison with quark mass functions measured 
in LQCD (see, e.g., Ref. |28|) can be made, one has also 
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to introduce a nonlocal dynamical modification of the 
Dirac vector part of the quark propagator for describ- 
ing the quark wave function renormalization as was re- 
cently accomplished within the nonlocal PNJL model in 
Refs. H|33|. 

Alternatively, one can start from QCD Dyson- 
Schwinger equations (DSEs), apply a symmetry- 
preserving truncation scheme, and solve the resulting 
equations for the Schwinger functions. In the recent 
years, this approach has reached a new level of matu- 
rity (see, e.g., Refs. (3ll - [33j for reviews). This ongoing 
progress in DSEs also includes ab initio type of calcu- 
lations [H|, H|| pertaining to PL and related to other 
functional approaches. 

There are attempts to derive nonlocal PNJL or PDSE 
models as a low-energy limit of QCD [36j], which used 
insights that have been obtained in the functional renor- 
malization group method [13, IH[. Functional methods 
also suggest close relation between the chiral and the con- 
finement transition [38l . |39| , and even the relation of both 
of them with the Ua(1) symmetry restoration [4(J. These 
studies are worth a systematic further development. 

Since QCD applications at T > are in general con- 
siderably more difficult than at T = 0, it is often useful 
to take a more phenomenological approach: one chooses 
a suitable (e.g., phenomenologically successful) effective 
dressed-gluon propagator model and solves the resulting 
T > DSEs for the Schwinger functions in the quark 
sector [32| . Such an approach has many advantages over 
NJL-type models as, e.g., a dynamical confinement mech- 
anism. However, as has been analyzed in [4l[ the criti- 
cal temperatures obtained within such non-perturbative 
low-energy QCD models turn out to be too low when 
compared with LQCD results. 

In the present work we suggest that this shortcoming 
might have its origin in residual color correlations and 
we investigate a generalization where the PL is coupled 
to the quark sector DSE (PDSE). For the purpose of 
this exploratory approach, we will restrict ourselves in 
this PDSE model to a rank-2 separable form of the effec- 
tive gluon propagator [42l | which at the rainbow-ladder 
level of truncation is equivalent to a full DSE model with 
translation-invariant gluon propagator once the model 
form factors of the interaction are chosen appropriately 
[43l | . For the time being one may take the pragmatic 
approach of modeling form factors of the rank-2 sepa- 
rable interaction such that both, the dynamical quark 
mass function and the wave function renormalization as 
measured in LQCD at zero temperature can be repro- 



duced to high accuracy. With an interaction model fixed 
in this way, the approach developed here can be used to 
predict thermodynamics and hadron properties of low- 
energy QCD at finite temperature, calibrate the results 
with LQCD and extend the approach subsequently to the 
whole QCD phase diagram, i.e., into regions presently in- 
accessible to LQCD. 

In the present work, by analyzing the temperature be- 
havior of the order parameters in the model, the dynami- 
cally generated light and strange quark masses as well as 
the PL variable, we find that the critical temperatures for 
the chiral and the deconfinement transitions measured by 
the peaks in the corresponding susceptibilities (defined 
here as the temperature derivatives of the order parame- 
ters) coincide at the per mille level of accuracy. We will 
discuss the effect of rescaling the critical temperature pa- 
rameter To of the PL potential p4l - l47j once applications 
with a finite number of quark flavors and a possible che- 
mical potential are considered and how this affects the 
width of the QCD transition region. 

As a possible application of this class of models we 
discuss the temperature dependence of scalar and pseu- 
doscalar meson properties at finite temperatures towards 
the chiral symmetry restoration. We devote special atten- 
tion to the investigation of a Gcll-Mann-Oakcs-Renncr 
(GMOR)- like relationship for the mass of pseudoscalar 
mesons at finite temperatures. We find that the GMOR- 
like relation proves to be robust up to the critical tem- 
perature, as a manifestation of the implemented confine- 
ment mechanism at the rainbow-ladder level of the de- 
scription, where the coupling to the PL provides a very ef- 
fective suppression of the quark degrees of freedom in the 
medium which otherwise induce a medium dependence 
which is unphysical below T c . In this regime, medium ef- 
fects have to be absent as long as hadronic excitations in 
the medium are disregarded. 



II. SEPARABLE PDSE MODEL 

A. Thermodynamical potential and order 
parameters 

The central quantity for the analysis of the thermody- 
namical behavior is the thermodynamical potential which 
in the PDSE approach is a straightforward generalization 
of the standard CJT functional [H HH 
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where the full quark propagator for the flavor / = u,d, s, 

= ij-pA f ((p%) 2 ,T) 

+i^ n Cf((p%) 2 ,T) + B f ((p%) 2 ,T)l2) 



is defined by the DSE for the quark self-energy E, see 
below. The Polyakov-loop potential is first taken in the 
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form [5C 
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with o(T) = a + ai(T /T) + a 2 {T /T) 2 , fc(T) - 
bs(To/T) . The corresponding parameters are taken 
from Ref. [50|, a = 3.51, oi = -2.47, a 2 = 15.22 and 
63 = —1.75 where they have been adjusted to fit the pres- 
sure obtained in lattice gauge theory simulations of SU(3) 
Yang-Mills theory. In most of the literature on the PNJL 
model, the parameter To = 270 MeV has been taken over 
for applications in QCD with Nf quark flavors while fol- 
lowing Ref. [45| its dependence on quark flavors and che- 
mical potential should be invoked. Accordingly, for the 
case Nf = 2+1 discussed in the present work, in [45[ 
the value To = 187 MeV is suggested with an error mar- 
gin of about 30 MeV. Applying the Matsubara formal- 
ism of finite temperature field theory, the squared quark 
4-momenta are to be replaced by (p™) 2 = (w") 2 + p 2 , 



lu" = Lu n + a(f>3, where uj„ = (2n + l)7rT are the fermionic 
Matsubara frequencies and the indices a = —1,0,4-1 
specify the three quark colors and their coupling to the 
parameter ^3 of the temporal gauge field. At vanishing 
chemical potential, the Polyakov loop is given by 
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In order to check the sensitivity to various paramcter- 
izations of the Polyakov-loop potential, we will also try 
the polynomial form 17 [ 
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with the temperature-dependent coefficient 
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b 2 {T) = a + oi 
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and the the set of parameters from Ref. [Uj, ao = 6.75, 
01 = -1.95, a 2 = 2.625, a 3 = -7.44, 63 = 0.75, and 
b 4 = 7.5. 

For the effective gluon propagator in a Feynman-likc 
gauge, g 2 Df u {p - q) = 5 m D(p 2 1 q 2 ,p ■ q), we employ a 
rank-2 separable ansatz [42( 

D(p 2 ,q 2 ,p-q) = DMp 2 )Mq 2 ) 

+D 1 T 1 (p 2 )(p-q)F 1 (q 2 ) , (7) 

so that the quark propagator amplitudes are given by 

B/((p«) 2 ,T) = m f + b f (T)M(pZ) 2 ) , (8) 
Af((p a n f,T) = l + a f {T)Fi{{pl) 2 ) , (9) 
Cf((p a n ) 2 ,T) = l + CfpmW) > (10) 

and their analytic properties are defined by the choice 
of the form factors. In the present work we will use the 
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(11) 
(12) 



which satisfy the constraints Jyj(0) = J-i(0) = 1 and 
J-"o(cxd) = .Fi (00) = 0. Their functional form can be cho- 
sen such that the 4- momentum dependence of the dynam- 
ical mass function M(p) = B(p)/A(p) and the the wave 
function renormalization Z(p) = 1/A(p) is in good agree- 
ment_[29j with LQCD simulations of the quark propaga- 
tor [23|. Models which employ a rank-1 separable ansatz 
(see, e.g., [H, [27]]) result in A(p) = Z(p) = 1 and miss 
an important aspect of quark dynamics in QCD. 

The temperature-dependent gap functions a/(T), 
bf(T) and c/(T) are obtained as solutions of the DSE 
for the quark self energy in rainbow-ladder truncation as 

m 
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where the denominator function is d/((p") 2 ,T) = 
^4((j£) a ,T) + «) 2 C 2 ((K) 2 ,T) + 5 2 ((K) 2 ,n 
Eqs. (|13 p -(|15 p correspond to minima of the ther- 
modynamical potential (H|) with respect to a vari- 
ation of the temperature dependent gap functions 
a/(T), bf(T), Cf(T) and have to be supplemented by a 
corresponding gap equation for the Polyakov loop which 
follows from the extremum condition dtt/d$> | m j n = 0. 
Once the gap equations are solved for different tem- 
peratures, one can extract the pseudocritical tempera- 
tures for chiral and deconfinement transitions from the 
peak positions of the temperature derivatives of the cor- 
responding order parameters, the quark mass functions 
m/(T) = [m° f + b(T)]/[l + af(T)] and the Polyakov loop 



where the index / (or /') stands for the quark (or an- 
tiquark) flavor (u, d or s), P is the total 4-momcntum, 
and t± = I ± P/2. The meson mass is identified from 
A(P 2 = -A/ 2 ) = 1. 

For example, with the separable interaction, the al- 
lowed form of the solution of Eq. (fTB}) for the pseudoscalar 
Bethe-Salpeter amplitude is 

T P (£; P) = 75 (iE P {P 2 ) + fF P (P 2 )) M? 2 )- (17) 

For scalar mesons we will use a truncated form of the 
Bethe-Salpeter amplitude (i.e. we take only the dominant 
contribution) 

T s (t ] P)=E s (P 2 ) F Q (l 2 ). (18) 

At T ^ in the Matsubara formalism, the 0(4) 
symmetry is broken by the heat bath and we have 
P — > P m = (v m , P) where v m = Im-nT . Bound states and 
the poles they generate in propagators may be investi- 
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$(T), respectively. For a discussion of the quark mass 
function as an order parameter of the chiral transition 
see, e.g., Refs. [53T4551 ] . 

B. Pion, kaon and sigma meson at finite 
temperature 

At T = the mass-shell condition for a meson as a 
qq' bound state of the Bethe-Salpeter equation (BSE) is 
equivalent to the appearance of a pole in the qq' scatter- 
ing amplitude as a function of P 2 . The qq' meson Bcthc- 
Salpctcr bound-state vertex Tff> (p, P) is the solution of 
the BSE 



(16) 

I 

gated through polarization tensors, correlators, or Bcthc- 
Salpctcr eigenvalues. This pole structure is characterized 
by information at discrete points v m on the imaginary 
energy axis and at a continuum of 3-momenta. One may 
search for poles as a function of P 2 thus identifying the 
so-called spatial or screening masses for each Matsub- 
ara mode. These serve as one particular characteriza- 
tion of the propagator and the T > bound states. In 
the present context, the eigenvalues of the BSE become 
A(P 2 ) — > A(^,-P 2 ;T). The spatial screening masses are 
identified by zeros of 1 - A(0, P 2 ;T). 

The general form of the finitc-T pseudoscalar and 
scalar Bethe-Salpeter amplitude allowed by the separable 
model for the lowest Matsubara mode = (as required 
for the spatial meson modes of interest here) is 

MC-P) = 75 (iE P (P 2 )+i-PFp(P 2 )) Jo((C) 2 ) 

r s (CP) = Es(P 2 )M(0 2 ) (19) 

One can then write the BSE for the spatial masses as 
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For example, BSE for the scalar a meson is 
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where uaj = Af/df, oqj = Cf/df and <tb,/ = Bf /df. 
Further details on the analysis of mcsonic spatial screen- 
ing masses in other channels and for specific model in- 
teraction kernels can be found, e.g., in Ref. (42l. l5ll l52j. 
Here, we have generalized this approach by accounting 
for the PL phase factors entering the quark propagators 
and the interaction kernels. 



III. RESULTS AND DISCUSSION 

For the numerical calculations we fix the free parame- 
ters of the model at T = as in Refs. [5lll5^.[5a |. to repro- 
duce in particular the vacuum masses of the pscudoscalar 
and vector mesons, M n = 140 MeV, M K = 495 MeV, 
M p = 770 MeV, the pion decay constant f n = 92 MeV, 
and decay widths, T p o^. e + e - = 6.77 keV, Tp^,^ = 151 
MeV as basic requirements from low-energy QCD phe- 
nomenology. 

For clarity we point out that at T = 0, the model 
without PL coincides with the otherwise same model 
including PL (just as they do in the NJL vs. PNJL 
case). We thus obtain the same parameter set as in 
Refs. [H [H [13, namely ?ti° = m° d = m° = 5.49 MeV, 
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115 MeV, D A 2 



m d = m " 
219, L»iA* = 



69, Ar 



0.758 GeV, Ai = 0.961 GcV and p Q = 0.6 GeV. 



A. Order parameters for chiral and deconfinement 
transition 

In Fig. [1] we show the resulting temperature depen- 
dence of the derivatives of the quark mass being an or- 
der parameter of the chiral phase transition and of the 
Polyakov loop expectation value as an order parameter 
of the deconfinement transition. The peak values are at- 
tained at the corresponding pseudocritical temperatures 
for the chiral (T x ) and deconfinement (T^) transitions, 
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FIG. 1: (Color online) Quark mass susceptibilities (blue 
dashed line: light flavors; red dash-dotted line: strange flavor) 
with coupling to the Polyakov loop (left panel) and without 
it (right panel) as a function of the temperature. Note that 
without coupling to the Polyakov loop the chiral transition 
temperature is unrealistically low and the peak value for the 
light flavors is different from that for the strange one. 



respectively. In the right panel of Fig. [T]we show the re- 
sults when the quark and gluon sectors are uncoupled. In 
this case we have in the light quark sector T x = 128 MeV, 
whereas Td = 270 MeV according to the parametrization 
of the PL potential in the pure gauge sector. The value 
obtained for T x is in the typical range found in the DSE 
approach [4l| . The peak position of the chiral suscepti- 
bility in the strange quark sector does not coincide with 
the one in the light quark sector in this case. When the 
quark and gluon sectors are coupled these temperatures 
get synchronized so that T c = T x = T d = 195 MeV, as is 
demonstrated in the left panel of Fig. [TJ 

At the same time, when coupling the PL potential to 
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FIG. 2: (Color online) Same as Fig. [TJ but as a function of 
the scaled temperature T/T c with T c = 195 MeV (left panel) 
and T c = T x = 128 MeV (right panel). Without coupling to 
the Polyakov-loop Td — 2.11 T c is outside the range shown. 



the chiral quark sector, the width of the transition region 
collapses to a tiny temperature interval around T c , as 
is demonstrated in Fig. [5] where the susceptibilities are 
shown as functions of the scaled temperature T/T c in the 
same interval with (left panel) and without (right panel) 
coupling the quark sector to the Polyakov loop potential. 

Both effects of coupling the chiral quark sector to the 
PL, the synchronization of the chiral and deconfmement 
transitions as well as the narrowing of the width of the 
QCD transition region, are obtained in a similar way for 
the polynomial PL potential J5]). 

The value obtained for the QCD transition tempera- 
ture, T c = 195 MeV (193 MeV) for the logarithmic (poly- 
nomial) PL potential, is closer to recent LQCD results 
than the one obtained in PNJL or rank-1 separable non- 
local PNJL models but unsatisfactory for a quantitative 
description. Within the framework of the PQM model, it 
has been suggested [45| to rescale the To parameter of the 
PL potential depending on the quark flavor content of the 
system and the chemical potential. We will follow such a 
prescription also in the present approach. In Fig. [31 we 
show the resulting temperature dependence of the order 
parameters for chiral symmetry breaking (the normal- 
ized mass function ra(T)/m(0)) and for deconfmement 
(the PL <J?(T)) for three values of T . According to [H| 
the case T = 187 MeV corresponds to Nf = 2 + 1 while 
To = 270 MeV is the value for the pure gauge theory 
where the deconfmement is a first order phase transition. 
The coupling to the chiral quark dynamics changes the 
character of this transition to a crossover. Lowering the 
To parameter to 187 MeV changes both deconfmement 
and chiral restoration to strong first order phase transi- 
tions! This change of character happens at the critical 
value T = 210 MeV, also shown in Fig. H 

In Fig. m we summarize this finding by showing the 
dependence of T c on the Tq parameter of the PL po- 
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FIG. 3: (Color online) Temperature dependence of the order 
parameters for chiral symmetry breaking (m(T)/m(0), blue 
lines) and for deconfmement ($(T), black lines) for different 
choices for the parameter To in the Polyakov-loop potential. 
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FIG. 4: (Color online) Pseudocritical temperature for the 
chiral restoration transition vs. parameter To of the Polyakov- 
loop potential in the logarithmic form ([3| (black circles) and 
in the polynomial form (blue squares). For further details, 
see text. 



tcntials ([3J and j5]). For the logarithmic PL potential 
([3]), the positions of first order transitions are charac- 
terized by the full dots connected by a solid line, while 
the crossover transitions are given as open dots con- 
nected by a dashed line. Two regions of linear depen- 
dence can be identified when using the logarithmic PL 
potential: T c = const + 0.30 T for T < 210 MeV and 
T c = const + 0.40 T for T > 210 MeV. When using 
the polynomial form of the PL potential, we find the 
linear dependence as T c = const + 0.36 Tq. The change in 
the character of the QCD transition from a crossover for 



T > 210 MeV to a first oder transition for T < 210 MeV 
is accompanied by a sudden change in slope at T = 210 
MeV. It is remarkable that the To- rescaling introduced 
to account for a quark flavor dependence of the PL poten- 
tial when applied to the nonlocal separable PDSE model 
considered here, results in an obvious contradiction with 
LQCD concerning the character of the QCD transition: 
while in LQCD for Nf = 2 + 1 the finite-T transition 
is a crossover [HI, [5?| , the application of the suggested 
reparamctrization with the corresponding value To = 187 
MeV leads in the present model to a first order transition. 

On the other hand, for the polynomial PL potential 
([5]) the transition is a crossover for any of the considered 
values of To. Fig. [4] illustrates this by the dashed line 
connecting the points depicted by squares. 



B. Meson screening masses at finite T 

Following the approach to spatial meson screening 
masses developed in [42j in its generalization by the cou- 
pling to the Polyakov loop as given above, we have eval- 
uated the temperature dependence of scalar and pseu- 
doscalar meson masses. In Fig. we show the results 
with (lower panel) and without (upper panel) PL cou- 
pling together with the behavior of the threshold to the 
continuum estimated by the sum of the corresponding 
quark mass functions. Since the PL coupling leads to a 
strong suppression of thermal quark excitations below 
the critical temperature, the continuum thresholds are 
almost constant with a sudden drop in the vicinity of T c . 
This behavior is reflected in the temperature dependence 
of the meson screening masses. As a quantitative mea- 
sure for the width of the QCD transition, we suggest to 
consider either the difference of a- and w- mass squared, 
A^._ ff = — M%, or the difference between the tempe- 
rature To-_27r where the a meson mass equals the double 
pion mass (the threshold for closing the a — > 2ir decay 
channel and the temperature T a — 2q where the a me- 
son mass equals the double quark mass (the threshold for 
opening the a — > 2q) decay channel). Both measures re- 
veal that the width of the transition region reduces from 
about 10 % without to about 1 % with PL coupling. Note 
that our results for the temperature dependence of the 
a- and tc- meson without PL coupling are very similar to 
those obtained earlier within the DSE approach [59| . 

A crucial part of the correct chiral behavior of the light 
pseudoscalar meson octet in the DS approach, where the 
pseudoscalars are both qq bound states and (almost-) 
Goldstone bosons of the dynamically broken chiral sym- 
metry of QCD, is the linear dependence of the squared 
pseudoscalar masses on the current quark masses m q as 
the chiral limit is approached 
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FIG. 5: Temperature dependence of pseudoscalar (71", K) and 
scalar (a) meson masses in the present model without (upper 
panel) and with (lower panel) coupling to the Polyakov-loop 
potential. The temperature T CT _2ir (T CT — 2q) where the a me- 
son mass equals the double pion mass (double quark mass) 
determines the threshold for closing (opening) the a — > 2n 
(a — > 2q) decay channel. 



M]f, the variable 



the present approach in the vacuum and at finite tem- 
peratures. 

To this end, one calculates M P 
"pion" mass, for different values of current light quark 
mass, m q , while the current strange quark mass is kept 
fixed. In Fig. [S] we show that this relation is very well ful- 
filled in the vacuum (T = 0) up to current quark masses 
well exceeding 10 m°. At finite temperatures, the GMOR- 
like relation (|22|) qualitatively holds well up to T ~ T c . 

From the temperature independence of the pion mass 
up to T c , together with the validity of the GMOR-likc 
relation in this range of temperatures, we can conclude 
that the temperature dependence of the chiral condensate 
must be mirrored by that of f%. It is a question of utmost 
importance for the phenomenology of hadronic matter 
whether such a statement would also hold when one goes 
beyond the rainbow-ladder level of description to which 
we restricted ourselves in the present work. 



M P = M) r 



const (m / + mfi ) 



(22) 



IV. CONCLUSIONS 



It is therefore interesting to investigate the validity of this 
Gell-Mann-Oakes-Renner (GMOR)-typc relationship in 



We have employed a Polyakov-loop Dyson- Schwinger 
equation approach to investigate the pseudocritical tern- 
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FIG. 6: Squared pseudo-scalar mass vs. current quark mass 
for different temperatures (data symbols, connected by lines 
to guide the eye), show a GMOR-like behavior which holds in 
a wide range of quark masses, and up to temperatures very 
close to the critical one. 



peratures for the chiral and deconfinement transitions 
for Nf = 2 + 1 quark flavors using a rank- 2 separable 
model for the effective gluon propagator. We find that 
the pseudocritical temperature T x = 128 MeV for the 
chiral restoration and that for deconfinement, Td = 270 
MeV, differ by more than a factor of two when the quark 
and gluon sectors are considered separately. But when 
the coupling is switched on these transitions get syn- 
chronized and the pseudocritical temperatures become 
coincident T c = T x = T d = 195 MeV. 

We have investigated the dependence of T c on the pa- 
rameter To of the Polyakov-loop potential. For the log- 
arithmic potential (|3|), we found two regions of linear 
dependence with a change in slope at To = 210 MeV, 
accompanied with a change of the character of the QCD 
transition from a crossover for To > 210 MeV to a first 
oder transition for T < 210 MeV. It is a remarkable find- 
ing of the present work that the To- rescaling to account 
for a quark flavor and chemical potential dependence of 
the PL potential, which was suggested in [45| and inves- 
tigated in greater detail in [46| for different parameteri- 
zations of the PL potential, when applied to the nonlocal 
separable PDSE model considered here, results in an ob- 
vious contradiction with LQCD concerning the character 
of the QCD transition. While in LQCD for N f =2 + 1 
the finite-T transition is a crossover, in the present model 
with T) = 187 MeV it is a first order transition. Never- 
theless, for a different form of the PL potential, namely 
the polynomial form ([5|), we find that the QCD transition 
remains crossover even for the smallest considered values 
of T . 

As a consequence for possible phenomenological ap- 
plications of the presented approach we discussed that 
the coupling of the Polyakov-loop to the chiral quark dy- 
namics narrows the temperature region in which chiral 
symmetry is approached. Quantitative measures for this 



region are the a-ir squared mass difference (M% — M%) and 
the difference of temperatures for opening the a — > qq de- 
cay channel (T a -q q ) and the closing of the a — > 2-7T decay 

The narrowness of the QCD transition region (AT/T C ) 
obtained from these measures is at the one-percent level 
and thus much too small for an adequate description of 
the QCD transition as obtained in recent LQCD studies 
(see, e.g., Ref. [60|). 

We conclude that the separable PDSE approach pro- 
vides an essential improvement of the chiral quark dy- 
namics in PNJL models and nonlocal PNJL models which 
use a rank-1 separable ansatz for the quark interaction 
kernel, since it provides a running of both, the dynamical 
quark-mass function and the wave-function renormaliza- 
tion in close agreement with LQCD simulations of the 
quark propagator. It also provides the strong-coupling 
aspect of a dynamical confinement mechanism due to the 
absence of real quark mass poles. 

However, the investigation of the temperature depen- 
dence of the chiral and deconfinement order parameters 
characterizing the (pseudo-) critical temperature and the 
width of the QCD transition reveals also some inadequate 
aspects of the present level of description of this transi- 
tion. The critical temperature is too high and the tran- 
sition region is too narrow when compared with LQCD 
results. A rescaling of the PL-potential results in a lower 
value for T c , in accordance with recent LQCD results, 
but at the price of a narrowing of the QCD transition re- 
gion, for the logarithmic PL potential even changing the 
character of the transition to a first order one, in striking 
contradiction with LQCD. 

We expect that going beyond the rainbow-ladder level 
by including hadronic fluctuations beyond the mean field 
[47l l6ll . [r33| will entail an improvement of the approach. 
As has been demonstrated recently by including tt and a 
fluctuations in a consistent 1/N C scheme [lH, [63[, going 
beyond the mean field will lead to a lowering of the chi- 
ral transition temperature. The width of the transition 
region, however, appears as a sensitive constraint for the 
choice of an appropriate functional form of the PL poten- 
tial. Its possible dependence on the inclusion of hadronic 
correlations deserves a detailed study. We plan to extend 
our work in this direction. 
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